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We  study  the  number  of  pivots  required  by  the  primal  network  simplex 
algorithm  to  solve  the  minimum-cost  circulation  problem.  We  propose  a  pivot 
selection  rule  with  a. bound  on  the  number  of  pivots,  for  an  n- 

vertex  network.  This  is  the  first  known  subexponential  bound.  The  network  sim¬ 
plex  algorithm  with  this  rule  can  be  implemented. to  run  in  n^oexy^  +  oo) 
the  special  case  of  planar  graphs,  we  obtain  a  polynomial  bound  on  the  number 
of  pivots  and  the  running  time.  We  also  consider  the  relaxation  of  the  network 
simplex  algorithm  in  which  cost-increasing  pivots  are  allowed  as  well  as  cost- 
decreasing  ones.  For  this  algonthm  we  propose  a  pivot  selection  rule  with  a  C'’  < 
bound  of  0(nm  -  min  {log(/lC);  mlopi])  on  the  number  of  pivots,  for  a  network 
with  n  vertices,  m  arcs,  and  integer  arc  costs  bounded  in  magnitude  by  C.  The 
total  running  time  is  0(nm  logn  ■  min  (OognC),  m  logn)).  This  bound  is  com¬ 
petitive  with  those  of  previously  known  algorithms  for  the  minimum-cost  circu¬ 
lation  problem. 
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1.  Introduction 

The  minimum-cost  circulation  problem  is  a  central  problem  in  network  optimization 
[15,30,35,40].  Much  research  has  been  devoted  to  finding  efficient  algorithms  for  the  problem. 
Classical  algorithms  for  the  problem  include  the  cheapest  augmenting  path  algorithm  [10,25,26], 
the  out-of-  kilter  algorithm  [17],  the  cycle-canceling  algorithm  [29],  and  the  primal  network  sim¬ 
plex  algorithm  [12,13,31].  Although  some  of  these  methods  perform  quite  well  in  practice,  they 
all  run  in  time  exponential  in  the  problem  size  in  the  worst  case.  Edmonds  and  Karp  [14]  were 

^  Research  panially  tupponei  by  the  Natianal  Science  Founditicn,  Grail  No.  DCR-MOS961,  and  the  Office  of  Naval 
Researdi.  Coniract  No.  N00014-r7-K04«7. 
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the  first  to  devise  a  polynomial-time  algorithm  for  the  problem.  Their  method,  which  can  be 
viewed  [33]  as  combining  the  cheapest  augmenting  path  technique  with  capacity  scaling,  runs  in 
0(m  logU  (m  +  n  logn))  time  *  on  a  network  containing  n  vertices  and  m  arcs  and  having  integer 
arc  edacities  of  magnitude  at  most  U. 

The  work  of  Edmonds  and  Karp  left  open  the  question  of  whether  tlie  problem  has  a 
strongly  polynomial  algorithm,  i.e.,  an  algorithm  with  a  running  time  bound  polynomial  in  the 
grajrfi  size  (m)  assuming  unit  time  per  ariLhmetic  operation,  as  well  as  polynomial  in  the  input 
size  (the  number  of  bits  needed  to  represent  the  graph  and  the  arc  capacities  and  costs)  assuming 
time  per  arithmetic  operation  polynomial  in  the  number  of  bits  in  the  operands.  Tardos  [39]  was 
the  first  to  discover  a  strongly  polynomial  algorithm.  Her  work  triggered  the  discovery  of  other 
polynomial  and  strongly  polynomial  algorithms  by  Orlin  [32,33],  Fujishige  [16],  Bland  and  Jen¬ 
sen  [6],  Gain  and  Tardos  [18],  Goldberg  and  Taijan  [20,21],  and  Ahuja,  Goldberg,  Orlin,  and 
Taijan  [2].  To  date,  the  best  strongly  polynomial  bound  known  is  Orlin’s  bound  of 
0(m  log/j  (m  +  «  logn))  [33].  For  networks  whose  arc  costs  are  integers  of  magnitude  at  most  C, 
Goldberg  and  Taijan’s  bound  of  0(.nm  log  (n^lm)  log(nC))  [20]  is  best  for  certain  ranges  of  n,m, 
and  C.  For  networks  whose  arc  capacities  and  costs  are  integral,  Ahuja,  Goldberg,  Orlin,  and 
Taijan’s  bound  of  Oinm  loglog  U  log  (nC))  [2]  is  best  for  certain  ranges  of  n,m,U,  and  C. 

These  polynomial  algorithms  depend  crucially  on  the  idea  of  scaling,  and  their  development 
left  open  the  question  of  whether  any  of  the  classical  algorithms  can  be  made  polynomial,  if  pos- 

AU  time  boundc  ve  tuted  for  •  unit-cost  random-access  machine  [I]. 
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sible  without  relying  on  scaling.  Oilin  [32]  provided  a  partial  answer  by  proposing  a  dual  net¬ 
work  simplex  algorithm  with  a  strongly  polynomial  time  bound.  His  result  does  use  scaling. 
Goldberg  and  Taijan  [21]  provided  another  answer  by  showing  that  Klein’s  cycle-canceling  algo¬ 
rithm  [29]  is  strongly  polynomial  if  a  mimmum  mean  cost  criterion  is  used  to  select  cycles  to 
cancel.  They  described  an  implementation  of  this  algorithm  with  a  running  time  of 
0(nm  logn  •  min  (log(rtC),  m  logn)).  This  algorithm  maintains  primal  feasibility,  although  the 
solutions  it  produces  are  not  necessarily  basic. 

A  remaining  open  question  is  whether  the  primal  rtetwork  simplex  algorithm  has  a  pivot 
selection  rule  that  guarantees  a  polynomial  or  strongly  polynomial  bound  on  the  number  of 
pivots.  As  Oriin  [32]  observed,  this  question  is  motivated  both  by  the  practical  efficiency  of  the 
network  simplex  algorithm  [3,23]  and  by  the  polyrwmial  bounds  that  have  been  obtained  for  the 
average-case  behavior  of  the  simplex  algoritiun  for  general  linear  programming  [9,38].  Recent 
work  of  Goldfarb  and  Hao  [22]  adds  immediacy  to  this  question.  They  considered  the  primal  net¬ 
work  simplex  algorithm  for  the  maximum  flow  problem  and  discovered  a  pivot  rule  guaranteeing 
0{nm)  pivots  and  O(n^m)  running  time.  With  tire  use  of  dynamic  trees  [36,37,40],  the  running 
time  of  their  algorithm  can  be  improved  to  0(nm  logn)  [19].  This  bound  is  within  a  factor  of 
logn  or  less  of  all  other  Imown  bounds  for  the  problem. 

In  this  paper  we  study  the  number  of  pivots  required  to  solve  the  minimum-cost  circulation 
problem.  We  propose  a  pivot  rule  that  has  an  bound  *  on  the  number  of  pivots  and 

the  same  bound  on  running  time  (with  a  different  additive  constant  in  the  expronent).  Although 

*  All  logarithnu  in  this  paper  are  base  two. 


this  bound  is  not  polynomial,  it  is  the  first  subexponential  bound  known.  For  the  special  case  of 
planar  graphs,  we  obtain  polynomial  bounds  on  the  number  of  pivots  and  the  running  time. 

We  also  consider  a  relaxation  of  the  primal  network  simplex  algorithm  in  which  cost- 
increasing  as  well  as  cost-decreasing  pivots  are  allowed.  For  this  algorithm,  we  propose  a  pivot 
rule  with  bounds  of  0(ww  •  min  {log(nC),  m  logn})  on  the  number  of  pivot::  and 
0{nm  logn  •  min  {log(nC),  m  logn})  on  the  total  running  time.  This  time  bound  is  the  same  as 
that  of  the  Goldberg-Taijan  cycle-canceling  algorithm.  This  is  no  coincidence;  our  results  are 
based  on  theirs.  In  particular,  the  notions  of  e-optimality  and  minimum  cycle  means  are  central 
to  our  work. 

This  paper  consists  of  five  sections  in  addition  to  the  introduction.  Section  2  contains  our 
network  terminology.  Section  3  presents  the  primal  network  simplex  algorithm  Section  4 
discusses  our  subexponential  bound  for  general  graphs  and  our  polynomial  bound  for  planar 
graphs.  Section  5  discusses  our  polynomial  bound  for  the  relaxed  algorithm.  Results  very  simi¬ 
lar  to  those  in  Section  5  have  been  obtained  independently  by  Orlin  (private  communication, 
1988)  and  Goldfarb  and  Hao  (private  communication,  1988)  using  essentially  the  same  methods. 
Section  6  contains  some  concluding  remarks. 

2.  Terminology 

Our  framework  for  studying  the  minimum-cost  circulation  problem  is  essentially  that  of 
Goldberg  and  Taijan  [20,21].  Let  G=(y,E)  be  a  directed  graph  with  vertex  set  V  of  size  n  ^2 
and  arc  set  £  of  size  m.  We  require  G  to  be  symmetric,  i.e.,  (v,^)  e  £  if  and  only  if  (h',v)  e  £. 
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Symmetry  implies  that  we  can  think  of  G  as  being  undirected  when  it  is  convenient  to  do  so.  We 
assume  without  loss  of  generality  diat  G  is  connected,  i.e.,  every  vertex  is  reachable  from  every 
other  vertex.  A  pair  of  arcs  (v.w),  (w.v)  is  called  an  edge  and  is  represented  by  the  unordered 
pair  {v.w}.  For  a  vertex  v,  we  denote  by  £(v)  the  set  {w  |  (v.h')  e  £}.  Graph  G  is  a  circulation 
network  if  each  arc  (v.w)  has  a  real-valued  capacity  u(v,w)  and  a  real-valued  cost  c(v,w).  We 
require  the  cost  function  to  be  antisymmetric,  i.e.,  c(v.h')  =  -c(>v,  v)  for  all  (v.w)  e  £.  We  denote 
by  C  the  maximum  absolute  value  of  an  arc  cost  if  aU  arc  costs  are  integers;  if  some  arc  cost  is 
nonintegral,  C  =  <». 

A  circulation  is  a  real-valued  function  / on  arcs  satisfying  the  foUowing  constraints: 

V  (v.w)  e  £,  /(v,w)  S  u(v,w)  (capacity  constraints):  (1) 

V  (v.w)  e  £,  /(v.w)  =  -/(w,v)  (flow  antisymmetry  constraints);  (2) 

V  V  €  V,  5^  /(v.w)  =  0  (conservation  constraints).  (3) 

w  c  £(v) 

The  cost  of  a  circulation  /  is  defined  by  the  following  formula: 

cost  if)  =  £  c(v,w)/(v.w)  (4) 

(v,w)e  E 
f(v,w)  >  0 

The  minimum-cost  circulation  problem  is  that  of  finding  a  circulation  of  minimum  cost. 

For  a  circulation  /  and  an  arc  (v.w).  the  residual  capacity  of  (v.w)  is 
Uf  (v.w)  =  «(v,w)  -/(v.w).  Arc  (v.w)  is  residual  if  M/(v,w)  >  0  arxl  saturated  if  u/  (v.w)  =  0.  We 
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denote  by  Ef  the  set  of  residua]  arcs  and  by  Gf  the  residual  graph  (V.E/).  An  edge  {v.h')  is  resi¬ 
dual  if  both  (v.w)  and  (h'.v)  are  residual;  otherwise  (v.h'}  is  saturated. 

A  cycle  is  a  sequence  (vo.vj),  (vj.vj).  Cv2.v3),...,(v*_i,vt)  of  arcs,  with  no  arc  repeated, 
such  that  Vo  =  V*.  The  cycle  is  simple  if  v;  ^  for  l<i  <j  <k.  The  residual  capacity  of  a  cycle 
is  the  minimuin  of  the  residual  capacities  of  its  arcs.  A  residual  cycle  is  a  cycle  whose  residual 
capacity  is  positive.  The  cost  of  a  cycle  is  the  sum  of  the  costs  of  its  arcs.  A  negative  cycle  is  a 
cycle  whose  cost  is  negative.  The  mean  cost  of  a  cycle  is  its  cost  divided  by  the  number  of  arcs  it 
contains.  The  minimum  cycle  mean  of  a  graph  with  arc  costs  is  the  minimum  of  the  mean  costs  of 
its  cycles,  and  a  minimutn  mean  cycle  in  such  a  graph  is  a  cycle  whose  mean  cost  is  minimum. 

The  following  well-known  theorem  characterizes  minimum-cost  circulations: 

Theorem  2.f  [11].  A  circulation  is  minimum-cost  if  and  only  if  there  are  no  negative  residual 
cycles. 

An  altcmaiive  characterization  of  minimum-cost  circulations  is  provided  by  linear  program¬ 
ming  duality  theory.  A  price  function  p  is  any  real-valued  function  on  the  vertices.  Given  a  price 
function  p,  the  reduced  cost  of  an  arc  (v.w)  is  Cp(y,w)  =  c(y,w)  +  p(y) -piw).  The  cost  of  any 
cycle  is  unaffected  by  such  a  transformation  of  the  arc  costs. 

Theorem  2.2  [15].  A  circulation /is  minimum-cost  if  and  only  if  there  is  a  price  function  p  such 
that 


V  (v.w)  €£,  m/v.w)  >  0  implies  Cpiv,w)  i  0  (optimality  constraints). 


(5) 
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We  need  a  notion  of  approximate  optimality  intnxluced  by  Tardos  [39]  and  independently 
by  Bertsekas  [43]-  This  notion  is  die  basis  of  several  minimum -cost  circulation  algorithms 
[2,5,20^1.39].  For  a  real  number  £^0,  a  circulation /is  t-optimal  if  there  is  a  price  fuiKtion  p 
such  that 

V  (v.w)  e  £,  Ufiv,w)  >  0  implies  Cp(v,w)  ^  -€  (e-optimality  constraints).  (6) 

Note  that  0-optimality  is  the  samff  as  optimality.  As  Bertsekas  discovered,  if  all  arc  costs 
are  integers,  then  e-optimality  for  a  sufficiently  small  e  is  enough  to  give  optimality; 

Theorem  23  [5]:  If  all  arc  costs  are  integers  and  e  <  1/n,  then  any  e-optimal  circulation  is 
optimal. 

The  key  to  our  results  is  the  connection  between  minimum  cycle  means  and  e-optimality, 
which  was  discovered  by  Goldberg  and  Taijan.  For  a  drculatian  /,  we  denote  by  p.(/)  the 
minimum  mean  cost  of  a  residual  cycle,  ix^  thetoinimum  cycle  mean  of  the  residual  graph  G/, 
and  by  e(/)  the  minimum  e  such  that  /is  e-optimaL 

Theorem  2.4  [20,21].  For  any  circulation/,  e(/)= max  [0, 

The  problem  of  computing  €(f)  is  dus  equivalent  to  die  problem  of  computing  the 
minimum  cycle  mean  of  Gf.  There  are  three  known  efficiect  algorithms  for  computing  minimum 
cycle  means:  Karp’s  0(nm )-time  algorithm  [27],  Karp  and  Orlin’s  0(nm  log«)-time  algorithm 
[28],  and  Orlin  and  Ahuja’s  log(«C))-time  algorithm  [34];  the  last  of  these  requires 

integer  arc  costs.  Given  an  e  such  that  /  is  e-optimal,  a  price  function  p  satisfying  the  e- 
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optiniality  constraints  can  be  found  in  0(nm)  time  using  the  Bellman-Ford  shortest  path  ligo- 
rithm,  as  discussed  by  Goldberg  and  Taijan  [20]. 

We  need  one  more  concept  related  to  the  ner>"orlc  simplex  algorithm.  A  circulation  /  is 
basic  if  the  set  of  residual  edges  fcnns  a  forest  (a  set  of  vertex-disjoint  trees). 

3.  Cycle-Canceling  and  the  Primal  Network  Simplex  Algorithm 

Theorem  2.1  suggests  the  following  cycle-canceling  algorithm,  discovered  by  Klein  [29], 
for  finding  a  minimum-cost  circulatioa  Begin  with  any  circulation.  Gf  there  exists  a  circulation, 
one  can  be  found  using  any  maximum  flow  algorithm.)  Repeat  the  following  step  until  there  are 
no  negative  residual  cycles:  Find  a  negative  residual  cycle  F  and  cancel  it  by  increasing  the  flow 
on  each  arc  of  F  by  an  amount  equal  to  the  residual  capacity  of  F.  Although  this  algorithm  can 
run  for  an  exponential  number  of  steps  if  a  bad  choice  of  cycles  to  cancel  is  made,  Goldberg  and 
Taijan  [21]  showed  that  if  a  minimum  mean  cost  residual  cycle  is  canceled  at  each  step,  the 
number  of  steps  is  polynomial,  specifically  0{nm  min  (log(nC),mlogn}).  A  related  rule  fcr 
choosing  cycles  to  cancel  can  be  implemented  to  that  the  total  running  time  of  the  algorithm  is 
0(nm  logn  •  min  (log(nC),m  logn}). 

The  cycle-canceling  algorithm  maintains  primal  feasibility  in  the  linear  programming  sense, 
i.e.,  it  maintains  a  circulation  at  each  step,  but  it  does  not  maintain  a  basic  solution,  i.e.,  a  basic 
circulation.  The  primal  network  simplex  algorithm,  henceforth  just  called  the  network  simplex 
algorithm,  is  a  variant  of  the  cycle-canceling  algorithm  that  does  maintain  a  basic  solution.  Tl.e 
simplex  algorithm  maintains  a  basic  circulation  / and  a  spanning  tree  T  whose  edges  are  edges  of 
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G,  such  that  T  contains  all  the  residual  edges.  An  arc  (v,w)  is  said  to  be  a  nee  arc  if  the 
corresponding  edge  (v.k')  is  in  T  and  a  nontree  arc  otherwise.  A  nontree  arc  (v.w)  defines  a 
basic  cycle  Y  (v.w)  consisting  of  the  arc  (v.w)  and  the  simple  path  of  tree  arcs  from  w  to  v.  Arc 
(v.w)  is  pivotable  if  Y  (v,w)  is  negative.  The  simplex  algorithm  consists  of  repeating  the  follow¬ 
ing  pivot  step  until  no  nontree  arc  is  pivotable;  Select  a  pivotable  nontree  arc  (v.w).  Cancel 
r  (v,w)  by  increasing  the  flow  on  every  arc  of  F  (v.w)  by  an  amount  equal  to  the  residual  capa¬ 
city  of  r  (v.w).  (This  amount  may  be  zero;  in  this  case  the  pivot  is  called  degenerate.)  Add  edge 
{v.w}  to  T  and  delete  some  edge  (x,y)  such  that  arc  (x,y)  is  on  F  and  (x,y)  is  saturated.  Arc 
(v.w)  is  called  the  eruering  arc  of  the  pivot  and  arc  (x.y)  is  called  the  leaving  arc,  it  is  possible 
for  the  two  to  be  equal.  The  pivot  is  said  to  be  a  pivot  on  (v.w). 

The  simplex  algorithm  differs  from  the  cycle -canceling  algorithm  in  that,  although  both 
algorithms  work  by  canceling  cycles,  the  simplex  algorithm  cancels  only  basic  negative  cycles, 
whereas  the  cycle-canceling  algorithm  cancels  ortly  residual  negative  cycles.  A  degenerate  pivot, 
i.e.,  one  canceling  a  nonresidual  basic  cycle,  does  not  change  the  circulation  /,  only  the  spanning 
tree  T.  The  relaxed  simplex  algorithm,  which  we  study  in  Section  5,  is  the  relaxation  of  the  sim¬ 
plex  algorithm  in  which  a  pivot  is  allowed  on  any  nontree  residual  arc  (v.w),  whether  or  not 
F  (v,w)  is  negative. 

Staitdard  implementations  of  the  simplex  algorithm  maintain  a  price  function  p  such  that 
every  tree  arc  has  a  reduced  cost  of  zero.  Such  a  price  function  can  be  defined  by  rooting  «hp  tree 
r  at  an  arbitrary  vertex  r  and  defining  p(v)  for  any  vertex  v  to  be  the  total  cost  of  the  tree  arcs 
along  the  simple  path  in  the  tree  from  r  to  v.  With  such  a  price  function,  a  nontree  residual  arc 
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(v.w)  is  pivotable  if  and  only  if  it  has  negative  reduced  cost.  Furthermore,  on  termination  of  the 
algorithm  the  circulation  / and  the  price  function  p  satisfy  the  optimality  constraints. 

The  simplex  algorithm  requires  an  initial  basic  circulation  to  get  started.  Such  a  basic  cir¬ 
culation  can  be  found  by  using  any  maximum  flow  algorithm  to  find  some  cinmlation  and  then 
repeatedly  finding  an  undirected  simple  cycle  of  residual  edges  and  increasing  (or  decreasing)  the 
flow  around  the  cycle  so  that  one  of  its  edges  becomes  saturated.  All  such  undirected  residual 
cycles  can  be  canceled  in  0(m  logn)  time  by  using  a  slight  variant  of  the  algorithm  of  Sleator  and 
Tarjan  [36]  for  making  a  flow  acyclic.  Having  found  a  basic  circulation,  a  suitable  initial  span¬ 
ning  tree  can  be  found  by  choosing  any  spanning  tree  containing  all  the  residual  edges.  Such  a 
tree  can  be  constructed  in  0(m)  time  using  graph  search. 

Remark.  The  idea  of  canceling  undirected  residual  cycles  can  also  be  used  to  convert  any 
minimum-cost  circulation  into  a  basic  minimum-cost  circulation.  Given  a  minimum-cost  circula¬ 
tion  /,  we  first  compute  a  price  function  p  with  respect  to  which  f  satisfies  the  optimality  con¬ 
straints.  This  t4Cves  0(nm)  time  as  noted  in  Seaion  2.  Any  residual  edge  |v,h')  is  such  that  both 
(v.H')  and  (w.v)  have  reduced  cost  zero.  We  then  cancel  undirected  residual  cycles  as  described 
above.  Canceling  any  such  cycle  does  not  change  the  cost  of  the  circulation,  since  any  such  cycle 
has  zero  reduced  cost.  All  the  cycle-canceling  takes  0(m  logn)  time.  □ 

The  possibility  of  degeneracy,  as  indicated  by  the  presence  of  one  or  more  saturated  edges 
in  the  spanning  tree  T,  requires  us  to  adjust  some  of  the  definitions  of  Section  2.  We  say  that  an 
arc  (v.w)  is  pseudo-residual  if  (v,^)  is  residual  or  (v.h')  is  in  T.  We  denote  by  E}  the  set  of 
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pseudo-residual  arc:  A  pseudo-residual  cycle  is  a  simple  cycle  consisting  of  pseudo-residual 
arcs. 

4.  A  Subexponential-Time  I^vot  Rule 

Our  pivot  rule  for  the  network  simplex  algorithm  is  based  on  the  cycle-canceling  algorithm 
of  Goldberg  and  Taijan.  Consider  th^  refinement  of  the  simplex  algorithm  that  consists  of  repeat¬ 
ing  the  following  siep  urul  there  is  no  negative  pseudo-residual  cycle;  Choose  a  minimum  mean 
cost  pseudo-residual  cycle  F  and  block  it  by  repeatedly  doing  pivots,  each  of  which  is  on  an  arc  of 
r,  until  some  arc  of  F  is  not  pseudo-refidual,  i.e.,  some  arc  of  F  is  a  saturated  nontrce  arc. 

We  postpone  a  discussion  of  how  to  choose  pivots  to  block  a  cycle  F.  First,  we  establish  a 
polynomial  bound  on  the  number  of  cycle  blockings.  Our  bound  is 
Oinm  ■  min  {log(nC),  m  logn)),  whi  h  is  the  same  as  the  bound  on  cycle  cancellations  for  the 
Goldoerg-Tarjan  algorithm.  The  derivation  of  our  bound  is  analogous  to  the  derivation  of  the 
Goldberg-Taijan  bound.  In  the  discussion  to  follow, /is  a  basic  circulation  and  T  is  a  spanning 
tree  containing  all  the  residual  edges.  The  analysis  requires  the  use  of  a  price  function  p  that  is 
ru)t  the  same  as  the  price  functiori  used  in  standard  implementations  of  the  network  simplex  algo¬ 
rithm.  (The  latter  price  fimetion  makes  the  reduced  cost  of  every  tree  arc  zero.)  This  distinction 
needs  to  be  remembered. 

For  a  real  number  e  i  0,  circulation  /  is  strongly  e-  optimal  if  there  is  some  price  function  p 
for  wnich  Cp(v,w)  S  -e  for  every  pseudo-residual  arc  (v.w).  Strong  e-optimality  implies  e- 
optimality  We  denote  by  £*(/)  the  minimum  e  for  which  /  is  strongly  e-optimal.  We  denote  by 


II* (/)  the  minimum  mean  cost  of  a  pseudo-residual  cycle.  Theorem  2.4  translates  into  the  follow¬ 
ing  result: 

Theorem4.1.  e*(/)  =  max  {0, 

Proof.  The  value  t(f)  is  a  function  only  of  the  set  of  residual  arcs  and  their  costs;  t(f)  is  the 
same  function  of  the  set  of  pseudo-residual  arcs  and  their  costs.  The  analogous  statement  is  true 
of  \L(f)  and  \L  (f).  The  theorem  is  thus  immediate  from  Theorem  2.4.  □ 

Let  t  =  t  (f),  and  let  p  be  a  price  function  with  respect  to  which  f  is  strongly  e-optimal. 
Holding  e  and  p  fixed,  we  study  the  effect  of  cycle  blockings  on  e*(/). 

Lemma  4 2.  Blocking  a  cycle  cannot  increase  e*(/). 

Proof.  (See  the  proof  of  Lemma  3.5  of  [21].)  Let  F  be  the  cycle  blocked.  Every  arc  of  F  has 
reduced  cost  exactly  -e.  Every  new  pseudo-residual  arc  (v.w)  created  by  a  pivot  step  during  the 
blocking  of  F  has  -eSCp(v,w)^e.  This  follows  by  cost  antisymmetry  and  the  fact  that,  just 
before  the  blocking,  every  tree  arc  (x.y)  satisfies  -c  S  Cp{x,y)  ^  e  by  strong  e-optimality.  Thus, 
after  the  blocking  of  F,  the  new  circulation  / is  still  strongly  e-optimal  with  respect  to  p  and  the 
new  T.  Therefore  e*(/)  S  e  after  the  blocking.  □ 

Lemma  43.  A  sequence  of  ml2-2n  -»■  2  Sm  cycle  blockings  reduces  e(f)  to  at  most  (l-l/n)  e. 

Proof.  (See  the  proof  of  Lemma  3.6  of  (21) .)  Consider  a  sequence  of  ml2-2n  +  2  cycle  block¬ 
ings.  Let  F  be  a  minimum  mean  cost  pseudo-residual  cycle  just  after  one  of  these  blockings. 
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(Cycle  r  might  or  might  not  be  the  next  cycle  to  be  blocked.)  Suj^se  F  contains  at  least  one  arc 
of  nonnegative  reduced  cost.  Let  /  ^  n  be  the  number  of  arcs  on  F.  By  strong  e-optimality, 
which  is  maintained  by  Lemma  4.2,  every  arc  of  F  has  reduced  cost  at  least  -e.  Thus  the  mean 
costof  Fisatleast-<(/-l)//)e  -(l-l//j)e.  This  implies  £*(/):$  (l-l/n)e. 

It  follows  that  if  any  of  the  ml2-2n  +2  cycles  blocked  contains  an  arc  of  nonnegative 
reduced  cost,  then  the  lemma  is  true.  (Once  e*(/)S(l-l/n)£,  this  inequality  remains  true  by 
Lemma  4.2.)  Suppose  that  each  of  the  m/2-2n  +  2  cycles  blocked  contains  only  arcs  of  negative 
reduced  cost.  Then  each  pivot  during  the  blockings  is  on  an  arc  of  negative  reduced  cost.  Furth¬ 
ermore,  each  cycle  blocking  produces  a  saturated  nontree  arc  (v.w)  of  negative  reduced  cost.  No 
later  pivot  can  be  on  either  (v,^)  or  (w.v).  Since  the  tree  T  always  contains  exactly  n-l  edges, 
there  are  always  exactly  2fl-2  tree  arcs.  Thus,  after  all  ot/2  -2«  +2  cycle  blockings,  every  non¬ 
tree  residual  arc  must  have  nonnegative  reduced  cost.  This  means  that  any  minimum  mean  cost 
pseudo-residual  cycle  contains  an  arc  of  normegative  reduced  cost.  By  the  argument  in  the 
preceding  paragraph,  t(f)£  (1-1  ln)t.  □ 

Theorem  4.4.  If  all  arc  costs  are  integers,  then  the  cycle  blocking  algorithm  terminates  after 
0(/wi  log(nC))  cycle  blockings. 

Proof.  (Sec  the  proof  of  Theorem  3.7  of  [21].)  Let /be  the  circulation  maintained  by  the  algo¬ 
rithm.  Initially  e*(f)  ^  C.  If  t(J)  <  1/n,  then  e(/)  =  0  by  Theorem  2.3,  since  strong  e-optimality 
implies  e-optimality.  By  Lemma  4.3,  each  group  of  m  cycle  blockings  reduces  e*(/)  by  a  faaor 
of  at  least  1-1 /n.  It  follows  that  the  algorithm  terminates  after  0(«m  log(nC))  cycle  blockings. 
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(Each  group  of  nm  cycle  blockings  reduces  e*(f)  by  at  least  a  constant  factor.)  □ 

To  obtain  a  strongly  polyrwmial  bound  on  the  number  of  cycle  blockings,  we  need  two 
more  definitions.  We  say  that  an  arc  (v.w)  is  t-fixed  if  f(y,w)  is  the  same  for  every  e-optimal  cir¬ 
culation  /.  The  following  theorem  of  Goldberg  arn!  Tatjan  generalizes  a  result  of  Tardos  [39]. 

Theorem  45  [20,21].  Let  e  >  0.  Suppose  a  circulation  is  e-optimal  with  respect  to  a  price  func¬ 
tion  p.  Suppose  further  that,  for  some  arc  (v.h'),  |  Cp(v.w)  j  ^ne.  Then  (v,w)  is  e-fixed. 

We  say  an  arc  (v,w)  is  dead  if  it  is  a  saturated  nontree  arc  and  its  flow  does  not  change  dur¬ 
ing  the  remainder  of  the  algorithm;  we  say  (v.w)  is  live  otherwise.  Every  arc  on  a  cycle  to  be 
blocked  is  live  since  each  such  arc  is  either  residual  or  a  tree  arc.  Thus  every  pivot  is  on  a  live 
arc. 


Theorem  4.6.  For  arbitrary  real-valued  arc  costs,  the  cycle  blocking  algorithm  terminates  after 
0(nm^  logn)  cycle  blockings. 

Proof.  (See  the  proof  of  Theorem  3.9  in  [21].)  Divide  the  cycle  blockings  into  groups  of 
k  =  cmn  logn,  where  c  is  a  suitably  large  positive  constant,  to  be  chosen  later.  We  shall  show  that 
each  group  of  k  cycle  blockings  converts  at  least  one  arc  from  live  to  dead.  The  theorem  foUows. 

Consider  a  group  of  k  cycle  blockings.  Let  / and  Z'  be  the  circulations  before  and  after  the  k 
cycle  blockings,  and  let  £  =  e*(/),  t'  By  Lemma  4.3,  e'  ^E/(2n)  for  c  a  sufficiently 

large  constant.  Let  F  be  the  first  cycle  blocked.  The  mean  cost  of  F  is  -c.  Let  p'  be  a  price 
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function  with  respect  to  which  f  is  strongly  e' -optimal.  Some  arc  of  F,  say  (v.w),  has  a  reduced 
cost  (with  respect  to  p')  of  at  most  -cS-2rte'.  By  Theorem  4.5,  (v,w)  is  e-fixed  after  all  the 
cycle  blockings.  Since  the  circulation  remains  strongly  c-optimal  by  Lemma  4.2,  the  flow 
through  (v,w)  cannot  subsequently  change.  Since  f  is  strongly  e' -optimal  with  respect  to p',  and 
since  the  reduced  cost  of  (v,^)  with  respect  to  p'  is  strictly  less  than  -e',  (v,w)  is  not  pseudo¬ 
residual,  i.e.,  it  is  a  saturated  nontree  arc.  Thus  (v,^)  is  dead  after  the  cycle  blockings.  This 
proves  the  claim.  □ 

Having  shown  that  the  number  of  cycle  blockings  is  polynomial,  we  must  still  devise  a  way 
to  perform  cycle  blockings  in  a  “small”  number  of  pivots.  Let  F  be  an  arbitrary  negative 
pseudo-residual  cycle.  The  following  method  will  block  F  in  at  most  pivots. 

We  need  to  work  with  condensations  of  the  graph  G.  If  (v,^)  is  an  arc  of  G,  the  condensed 
graph  G(v,w)  is  formed  by  identifying  vertices  v  and  w.  The  operation  of  forming  GCv,^)  is 
called  condensing  iv,w).  Note  that  GCv,^)  will  contain  loops  (arcs  of  the  form  (x,x))  and  may 
contain  multiple  arcs  (two  or  more  arcs  of  the  form  (x,y)).  We  allow  such  loops  and  multiple 
arcs  to  remain;  we  call  G(v,w)  a  multigraph.  If  5  is  a  subset  of  the  arcs  of  G,  we  form  the  con¬ 
densed  graph  G(S)  by  condensing  each  of  the  arcs  in  5  in  turn.  Up  to  renaming  of  vertices,  the 
order  of  these  condensations  is  irrelevant. 

Let  7  be  the  current  spanning  tree.  Let  p  be  a  price  ftmction  such  that  every  tree  arc  has 
zero  cost.  In  what  follows,  “reduced  cost”  means  with  respect  to  p\  the  reduced  cost  of  an  arc  in 
a  condensed  graph  is  defined  to  be  the  reduced  cost  of  the  corresponding  arc  in  the  original  graph. 
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We  begin  the  blocking  of  cycle  F  by  defining  a  sequence  of  triples 
(Gi.Fi.Fi).  (G2.F2,r2).—  (PkfTk.Tk)-  For  1  Sk,  Gi  is  a  condensation  of  G  fonned  by  con¬ 
densing  some  of  the  arcs  of  T;  is  the  spanning  tree  of  G,  whose  arcs  correspond  to  those  of  T\ 
r,  is  a  negative  reduced-cost  simple  cycle  in  G,  whose  arcs  correspond  to  those  of  F.  The  first 
triple  is  (G 1  .Tj  ,Fi )  =  (G,7,F)-  For  i  >  1,  the  i**  triple  is  defined  if  the  i-V  triple  is  defined  and 
I  Fj_i-r,_i  I  ^  2;  here  Fi_i  and  r,_i  are  interpreted  as  sets  of  arcs,  and  absolute  value  signs  denote 
set  size.  If  1  Fi_i-7,_i  \  >  2,  G,  =  G,_i  (7i_i  -F,_i).  The  tree  7,  is  the  spanning  tree  of  G,  consist¬ 
ing  of  the  nonloop  arcs  that  correspond  to  arcs  of  7,_i .  The  cycle  F,  is  chosen  to  be  any  negative 
reduced-cost  simple  cycle  in  G,  ,  all  of  whose  arcs  correspond  to  arcs  of  r,_i . 

Lemma  4.7.  For  2^i£k,a  suitable  cycle  Fj  can  be  chosen. 

Proof.  By  induction  on  i.  There  is  a  cycle  FJ  in  G„  not  necessary  simple,  whose  arcs  correspond 
to  those  of  Fj_i.  Cycles  F,_i  and  F,'  have  the  same  reduced  cost.  By  the  induction  hypothesis, 
F,_i  has  negative  reduced  cost;  thus  so  does  Fj.  The  arcs  of  Fj  can  be  partitioned  into  simple 
cycles,  at  least  one  of  which  must  have  negative  reduced  cost  and  hence  can  be  chosen  as  F,-. 

Since  cycle  F,  for  1  S  *  is  simple,  it  contains  a  loop  only  if  the  loop  is  the  entire  cycle. 
Thus  Fj  contains  no  zero  reduced-cost  loops,  which  means  that  any  arc  (v,^)  e  F,  n  7,  has  v  ^  w. 

Suppose  we  perform  a  pivot  on  die  arc  (v,w)  in  F  correspondirig  to  the  unique  arc  in  Ft-7*. 
Since  F*  has  rtegative  reduced  cost,  so  does  (v,w),  and  we  can  indeed  pivot  on  (v,w).  Let  (x,y)  be 
the  leaving  arc  of  the  pivot.  Let  j  be  maximum  such  that  x  and  y  are  not  condensed  in  Gj.  Then 
the  pivot  in  G  that  replaces  (v,w}  by  {x,y}  in  7 corresponds  to  a  pivot  in  each  G,  for  1  ^  J  that 
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forms  a  new  tree  r-  by  replacing  the  edge  corresponding  to  {v.^]  by  the  edge  corresponding  to 
{x,y}.  This  observation  leads  to  the  following  pivot  rule  for  blocking  F.  Construct  a  sequence  of 
triples  (G 1  .Fj  .Fi),  (G2.F2.F2),...,(G*.r*,Ft)  as  ^cified  above.  Then  repeat  the  following  step 
until  some  arc  of  F  becomes  a  saturated  nontree  arc: 

Pivot  and  Restore  Triples.  Pivot  on  the  arc  (v.w)  in  F  corresponding  to  the  unique  arc  in  Fi-7*. 
Let  (x,y)  be  the  leaving  arc  of  the  pivot.  Let  j  be  maximum  such  that  x  and  y  are  not  condensed 
m  Gj.  Perform  the  corresponding  pivot  in  G,  for  1  <  j  <j  to  form  T,  frcm  F,.  Restore  the 
sequence  of  triples  by  beginning  with  (Gi.F'i.Fi).-  (G^.F^.F^)  and  extending  the  sequence 
using  the  inductive  definition  above.  (The  new  number  of  triples  k  may  differ  from  the  old 
number.) 

In  order  to  measure  the  progress  of  this  algorithm,  we  define  a,  =  1  F,-F,  |  for  1  ^  i  <  k.  We 
define  the  signature  of  a  sequence  of  k  triples  to  be  the  sequence  ai,a2,...,0(t.  We  order  signa¬ 
tures  lexicographically,  i.e.,  ai,a2,...,a*  <  Pi,P2 . P;  ’f  either  there  is  some  j  S  min  {k.i]  such 

that  ce,  =  Pi  for  1  ^ i  < y  and  <  pj, or /t  <  / and  a,  =  Pi  for  l^i^k. 

Lemma  4.8.  Each  pivot  step  decreases  the  signature  or  teiminates  the  cycle  blocking. 

Proof.  Consider  a  pivot  on  an  arc  (v,h').  Let  ai,a2,....at  and  Pi,P2 . P/  be  the  signatures 

before  and  after  the  pivot,  respectively.  Let  (x,y)  be  the  leaving  arc,  and  let  j  be  maximum  such 
that  X  and  y  are  not  condensed  in  Gj.  If  j  <k,  then  the  arc  in  Gy  corresponding  to  (x,y)  is  in 
Fy-F y,  since  it  is  ore  of  the  arcs  condensed  to  form  Gy  +  j .  Furthermore,  for  1  ^  i  <  j,  the  arc  in 
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Gi  corresponding  to  (x.y)  is  in  T,  n  T,.  It  follows  that  a,  =  P,  for  <j  and  (Xj  <  py. 

Suppose  on  the  other  hand,  that  j  =  k.  Then  the  basic  cycle  in  G*  on  which  the  pivot  occurs 
is  r*.  and  arc  (x.y)  must  lie  on  T.  The  pivot  makes  (x,y)  a  saturated  nontree  arc,  which  ter¬ 
minates  the  blocking.  □ 

Lemma  4.9.  For  \  ^i^k,  CLi  =  \  Ti—Ti  |  ^  n/2'“* . 

Proof.  For  \  <i  ^k,  let  t,  be  the  number  of  vertices  v  in  G,  such  that  v  is  incident  to  some  arc  of 
7,-r,.  We  claim  that  a,  ^  t,  or  a,  =  1.  To  verify  the  claim,  suppose  that  a,  ^  2.  The  cycle  F, 
consists  of  arcs  not  in  T,  alternating  with  paths  of  arcs  in  7,  .  (Some  of  these  paths  may  consist  of 
single  vertices  and  no  arcs.)  Some  vertex  along  each  path  in  F,  n  7,  must  be  incident  to  an  arc  of 
7,-F,,  since  there  are  a,  ^  2  such  paths,  which  are  connected  in  7,.  Hence  a,  ^  T,. 

Now  we  show  that  T,  S  x,_i  /2  for  2^i  ^k.  This  inequality  implies  x,  <  n/2‘“’  for  1  S  i  <  it 
by  induction  on  i,  since  Xj  S  n.  This  in  turn  gives  the  lemma.  To  show  that  x,  S  x,_i  /2,  we 
observe  that  every  vertex  v  in  G,_i  incident  to  an  arc  in  7,_i  -F,_i  is  condensed  into  another  such 
vertex  in  forming  G,,  namely  the  other  vertex  incident  to  the  arc.  Furthermore  if  v  is  a  venex  in 
Ci  incident  to  an  arc  in  Ti-F/,  then  the  corresponding  vertex  in  G.-i  is  incident  to  the 
corresponding  arc,  and  this  arc  is  in  7,_i  -r,_i  because  every  arc  in  F,-  corresponds  to  an  arc  in 
r,_i.  It  follows  that  each  vertex  counted  in  x,  is  formed  by  condensing  two  or  more  vertices 
counted  in  x,_i ,  which  implies  that  x,  S  Xj_i  12.  □ 

Theorem  4.10.  The  cycle-blocking  algorithm  blocks  F  in  at  most  pivots. 
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Proof.  By  Lemma  4.8,  each  pivot  either  terminates  the  algorithm  or  decreases  the  signature. 

Do*"! 

Lemma  4.9  implies  that  there  are  at  most  11^  S  distinct  possible  signatures.  The 

theorem  follows.  □ 

Combining  Theorems  4.4, 4.6,  and  4.10,  we  obtain  the  main  result  of  this  section: 

Theorem  4.11.  The  network  simplex  algorithm  with  pivot  selection  based  on  minimum  mean 
cost  cycle  Mocking  terminates  in  ■  min  {log(r»C),  mlogn))  pivots  and 

(5(„0og»V2+ 1^2  .  {log(rtO,mIog«))  running  time. 

Proof.  The  bound  on  pivots  is  immediate.  It  is  not  hard  to  implement  the  algorithm  so  that  the 
amortized  time  *  per  pivot  is  0(m);  we  leave  this  as  an  exercise.  □ 

We  ctmjecture  that  the  amortized  time  F>er  pivot  can  be  reduced  to  0{n)  or  less,  but  we  shall 
not  pursue  this  issue  here,  since  the  bound  on  pivots  is  so  large.  Instead,  we  consider  the  case  of 
planar  graphs  and  show  that  for  such  graphs  the  number  of  pivots  can  be  made  polynomial. 

Suppose  that  G  is  pianar.  Planarity  implies  that  m  =  0(n)  [7].  We  shall  show  that  a  judi¬ 
cious  choice  of  and  Fs  in  the  cycle-blocking  algorithm  guarantees  that  every  signature  con¬ 
tains  at  most  three  terms,  i.e.,  one  of  ai,a2.  and  Oj  is  1.  This  implies  that  pivots  suffice  to 
block  a  cycle,  since  every  term  in  a  signature  except  the  last  is  between  2  and  n,  inclusive,  and  the 
last  is  1. 

By  "■mortued  time"  we  mean  the  time  per  operation  averaged  over  a  worst  case  sequence  of  operations. 

See  [41]. 
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Form  G\  from  Gi  =G  by  deleting  all  arcs  except  those  in  Ti  uFj.  Regard  Gj  as  an 
undirected  graph  and  embed  it  in  the  plane.  Assume  without  loss  of  generality  that  G\  is  embed¬ 
ded  so  that  Fi  is  directed  clockwise  around  its  interior.  Consider  the  boundaries  of  the  faces 
inside  F j ,  oriented  clockwise  around  their  interiors.  These  boundaries  are  simple  cycles  that  par¬ 
tition  the  arcs  of  Fj  and  the  tree  arcs  inside  Fj.  The  sum  of  the  costs  of  these  cycles  is  exactly  the 
cost  of  Fj;  thus  at  least  one  such  cycle  must  have  negative  cost.  Let  r\  be  such  a  negative  cycle. 
Choose  F2  to  be  any  negative  simple  cycle  in  G2  whose  arcs  correspond  to  those  of  F'l . 

If  IF2  -  7'2  I  ^  2,  proceed  as  follows  to  choose  F3.  Form  G2  from  G\  by  condensing  the 
arcs  of  T y  lying  inside  Fj  and  deleting  all  arcs  except  those  corresponding  to  noncondensed  arcs 
of  F'l  u  Fi .  Let  72  be  the  cycle  in  G2  corresponding  to  F'l ,  and  let  G2  be  the  spanning  tree  of 
G2  corresponding  to  the  noncondensed  arcs  of  Fi.  Condensing  the  planar  embedding  of  C'l 
gives  a  planar  embedding  of  C2  (regarded  as  an  undirected  graph)  such  that  F2  is  the  boundary  of 
a  face,  with  F2  oriented  clockwise  around  its  interior.  Consider  the  boundaries  of  the  other  faces 
of  G2.  oriented  counterclockwise  around  their  interiors.  These  boundaries  are  simple  cycles  that 
partition  the  arcs  of  H  Because  the  embedding  is  planar  and  7^  is  a  spanning  tree,  each 
such  cycle  contains  exactly  one  arc  of  n  -^2-  Let  F3  be  such  a  cycle  that  contains  a  negative- 
cost  arc  corresponding  to  an  arc  of  F2.  Choose  Fj  to  be  the  simple  cycle  in  G3  corresponding  to 
Then  IF3  -Fjl  =  1. 


Theorem  4.12.  If  G  is  a  planar  graph,  and  F2  and  F3  are  chosen  in  the  cycle-blocking  algorithm 
as  described  above,  then  the  number  of  pivots  per  cycle  blocking  is  at  most  n},  the  total  number 
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of  pivots  needed  to  find  a  minimum-cost  circulation  is  0(n*  min  {  log(nC),  nlogn}),  and  the 
time  needed  to  find  a  minimum-cost  circulation  is  0(n*  •  min  {log(nC),  nlogn)). 

Proof.  The  discussion  above  and  Lemma  4.8  imply  that  the  number  of  pivots  per  cycle  blocking 
is  at  most  n^.  The  bound  on  total  pivots  follows  from  Theorems  4.4  and  4.6  and  m  =  0(n).  The 
time  per  pivot  is  0(n)  if  a  linear-time  planarity  algorithm  [8,24]  is  used  to  embed  G.  (ActuaUy,  it 
suffices  to  embed  G  once  and  modify  the  embedding  as  needed  to  choose  entering  arcs  for  pivot 
steps.)  □ 

5.  A  Polynomial-Time  Relaxed  Network  Simplex  Algorithm 

The  results  in  Section  4  relate  the  efficiency  of  the  network  simplex  algorithm  to  that  of 
cycle  blocking.  If  we  relax  the  network  simplex  algorithm  to  allow  cost-increasing  pivots  in 
addition  to  cost-decreasing  ones,  then  cycle  blocking  becomes  much  easier.  The  following  argu¬ 
ment  shows  that  at  most  n  pivots  suffice  to  block  a  cycle.  Let  F  be  the  cycle  to  be  blocked,  and 
let  V  be  an  arbitrary  vertex  on  F.  The  rule  for  pivoting  is  to  pivot  on  the  first  nontree  arc  after  v 
on  F.  If  this  nile  is  used,  then  after  k  pivots  cither  F  is  blocked  or  the  first  k  arcs  after  v  on  F  are 
tree  arcs.  Thus,  after  n  or  fewer  pivots,  F  must  be  blocked.  This  gives  the  following  result: 

Theorem  5.1.  With  the  pivot  rule  for  cycle  blocking  proposed  above,  the  relaxed  netwoik  sim¬ 
plex  algorithm  takes  0(/t^m  •  min  {log(nC),  mlogn))  pivots  and  0(n^m^  ■  min  {log(nC), 


mlogn))  total  time. 
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Proof.  The  bound  on  pivots  foUows  from  nteorems  4.4  and  4.6  and  the  discussion  above. 
The  time  for  blocking  a  cycle  is  0(nm)  to  find  a  suitable  (minimum  mean  cost)  cycle  to  block 
plus  Oim)  for  each  of  the  n  pivots,  for  a  total  of  Oinm).  The  total  time  bound  follows  from 
Theorems  4.4  and  4.6.  □ 

We  shall  reduce  the  bourxl  on  pivots  in  Theorem  5.1  by  a  factor  of  n  and  the  bound  on  time 
by  a  factor  of  nm/logn.  To  do  this  we  first  relax  the  rule  for  choosing  cycles  to  block.  We  follow 
the  approach  Goldberg  and  Taijan  [21]  used  to  develop  a  more  efficient  version  of  the  cycle¬ 
canceling  algorithm. 

Recall  that  an  arc  (v.w)  is  pseudo-residual  if  it  is  residual  or  a  tree  arc,  and  that  a  basic  cir¬ 
culation  /  is  strongly  t-optimal  if  there  is  a  price  function  p  for  which  Cp(»',H')  >-c  for  every 
pseudo-residual  arc  (v.w).  The  new  algorithm  explicitly  maintains  such  a  price  function  p.  On 
Section  4,  such  a  price  ftmction  was  used  only  to  analyze  the  algorithm,  not  to  define  iu)  For  a 
basic  circulation  /  and  a  price  function  p,  we  denote  by  t{f,p)  the  minimum  e  >  0  such  that 
Cp(v.w)  t  -£  for  every  pseudo-residual  arc  (v.w).  We  call  a  pseudo-residual  arc  (v.w)  admissible 
if  Cp(v.v,>)  <  0.  The  idea  of  the  algorithm  is  to  repeatedly  perform  pivot  steps  on  admissible  arcs 
until  there  is  no  cycle  of  admissible  arcs,  and  then  modify  p  to  reduce  t  {f,p).  To  be  precise, 
given  an  initial  basic  circulation  /  and  a  spanning  tree  T  containing  the  residual  edges,  the  algo¬ 
rithm  first  computes  a  {Mice  function  p  such  that  t*(f)  =  c*(f,p).  Then  the  algorithm  repeats  the 
following  two  stepK  until /is  optimal: 

Step  1  (block  admissible  cycles).  Repeatedly  perform  pivots  on  admissible  arcs  until  there  is  no 
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cycle  of  admissible  arcs. 

Step  2  (tighten  prices).  Modify  p  so  that  t(f,p)  decreases  to  at  most  (1-1 /n)  times  its  former 
value. 

A  proof  like  that  of  Lemma  4.2  shows  that  Step  1  cannot  increase  t{f,p).  We  postpone  a 
discussion  of  the  implementation  of  Step  1  in  favor  of  Step  2. 

Step  2  can  be  implemented  eilfier  using  a  minimum  cycle  mean  calculation,  which  takes 
0(nm)  time  (see  Section  2),  or  with  the  following  simple  two-step,  0(m)-time  computation. 

Step  2a  (compute  levels).  For  each  vertex  v,  compute  a  level  L(v),  defined  recursively  by 
L(v)  =  0  if  V  has  no  incoming  admissible  arcs,  L(y)  =  1  +  max  (L(«)  |  (u,v)  is  an  admissible  arc) 
otherwise. 

Step  2b  (compute  new  prices).  Let  t  =  t*(f,p).  For  each  vertex  v,  replace  p(v)  by 
p(v)-(e/n)L(v). 

Lemma  52.  Either  implementation  of  Step  2  is  correa. 

Proof.  Lemma  4.2  of  [21]  states  that  Steps  2a  and  2b  correctly  implement  Step  2.  The  imple¬ 
mentation  using  a  minimum  eyrie  mean  calculation  n*duces  t(f,p)  to  t  (f),  which  is  the 
minimum  possible  value  of  t(f,p).  Since  Steps  2a  and  2b  will  reduce  t*(f,p)  to  at  most 
(1  -1  /n)  times  its  former  value,  so  must  the  minimum  cycle  mean  calculation.  □ 
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Suppose  we  implement  Step  2  by  using  Steps  2a  and  2b,  except  that  every  n"*  iteration  of 
Step  2  we  use  a  minimum  cycle  mean  computation  instead.  Then  the  amortized  time  per  iteration 
of  Step  2  is  0(m).  The  number  of  iterations  of  Steps  1  and  2  is  0(n  log(nC)).  This  is  because 
t  (f,p)  <  C  initially,  the  algorithm  terminates  when  t  (f,p)  <  l/n,  and  every  n  iterations  decrease 
t  (f,p)  by  at  least  a  constant  factor.  Furthermore,  a  proof  like  that  of  Lemma  4.6  shows  that  the 
number  of  iterations  of  Steps  1  and  2  is  0(nm  log/i).  Thus  we  obtain  the  foDowing  result; 

Theorem  5.3.  The  relaxed  network  simplex  algorithm  terminates  in 

0(n  ■  min  {log^nCl.^logn})  iterations  of  Steps  1  and  2.  The  amortized  time  per  iteration  of  Step 
2  is  0(m). 

It  remains  to  implement  Step  1.  We  shall  describe  a  way  to  carry  out  Step  1  in  0(m) 
pivots.  Then  we  shall  describe  how  to  extend  the  dynamic  tree  data  structure  of  Sleator  and  Tar- 
jan  [36,37,40]  so  that  it  can  be  used  to  perform  each  pivot  step  in  OQogn)  amortized  time  and  all 
of  Step  1  in  0(m  logn)  time. 

Our  implementation  of  Step  1  maintains  a  partition  of  the  vertices  into  two  classes:  marked 
and  unmarked.  Initially  all  vertices  are  urunarked.  A  vertex  v  becomes  marked  when  the  algo¬ 
rithm  discovers  that  there  is  no  path  of  one  or  more  admissible  arcs  from  v  to  an  unmarked  ver¬ 
tex.  The  algorithm  maintains  a  root  vertex  r  for  the  spanning  tree  T,  such  that  the  following 
invariant  holds; 


(*)  If  (v.w)  is  an  admissible  tree  arc  such  that  w  is  unmarked,  then  either  w  is  the  parent  of  v  in 
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T  (i.e.,  w  lies  on  the  tree  path  from  v  to  r),  or  (v.w)  has  always  been  an  admissible  tree  arc 
with  V  the  parent  of  w.  (In  particular,  if  the  latter  case  holds  then  (v,w)  was  not  added  to  T 
by  a  pivot.) 

Each  pivot  is  on  an  arc  of  the  form  (r.v).  The  algorithm  consists  of  the  following  steps: 

Step  la  (initialize).  Unmaric  aU  vertices.  Selea  any  vertex  r  and  root  T  at  r. 

Step  lb  (pivot).  If  there  is  no  admissible  arc  of  the  form  (r.v)  with  v  unmarked,  mark  r  and  go  to 
Step  Ic.  Otherwise,  let  (r.v)  be  such  an  arc.  If  (r.v)  is  a  tree  arc,  root  T  at  v,  replace  r  by  v,  and 
repeat  Step  lb.  If,  on  the  other  hand,  (r,v)  is  a  nontree  arc,  pivot  on  (r,v).  Let  ix,y)  be  the  leav¬ 
ing  arc  of  the  pivot.  Root  the  tree  at  x,  replace  r  by  x,  ?Jid  go  to  Step  Ic. 

Step  Ic  (reroot).  If  every  vertex  is  marked,  stop.  Otherwise,  let  v  be  any  unmarked  venex.  Let  x 
be  the  unmarked  vertex  closest  to  r  along  the  tree  path  from  v  to  r.  If  x  ^  r,  root  the  tree  at  x  and 
replace  r  by  x.  Go  to  Step  lb. 

Lemma  5.4  Steps  la-lc  correctly  implement  Step  1  and  maintain  invariant  (*). 

Proof.  Pivoting  on  an  admissible  arc  cannot  create  any  new  admissible  arcs.  Consider  Step  lb. 
It  there  is  no  admissible  arc  (r.v),  then  marking  r  preserves  the  invariant  that  there  is  no  path  of 
admissible  arcs  from  a  marked  vertex  to  an  unmarked  vertex.  If  (r.v)  is  an  admissible  tree  arc, 
then  rerooting  the  tree  at  v  obviously  preserves  (*).  If  (r.v)  is  an  admissible  nontree  arc,  lerooting 
the  tree  after  a  pivot  on  (r.v)  also  preserves  (*).  since  the  entering  arc  (r.v)  has  v  the  parent  of  r 
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aficr  the  reiooting  (unless  (r.v)  is  the  leaving  arc),  and  the  rerooting  leaves  invariant  the  parent 
and  the  child  vertices  of  all  other  tree  arcs. 

Rerooting  the  tree  as  in  Step  (Ic)  also  preserves  (*),  since  each  tree  arc  (y,z)  along  the  tree 
path  from  v  to  r  has  z  marked,  and  the  reversal  (z,y)  of  such  an  arc  has  y  the  parent  of  z  after  the 
rerooting.  □ 

Lemma  5  J.  The  number  of  iterations  of  Steps  lb  and  Ic  is  0(m). 

Proof.  We  need  only  count  iterations  of  Step  lb,  since  every  iteration  of  Ic  follows  an  iteration 
of  lb.  Each  iteration  of  lb  either  marics  a  vertex,  which  can  happen  at  most  n  times,  or  creates  a 
new  admissible  tree  arc  of  the  form  (x,y)  with  y  unmarked  and  y  the  parent  of  x,  or  does  a  pivot. 
Let  (x,y)  be  the  leaving  arc  of  a  pivot.  If  (x,y)  is  not  admissible  ory  is  marked,  (x,y)  will  never 
later  be  the  entering  arc  of  a  pivot.  If  (x.y)  is  admissible  arxl  y  is  unmarked,  then  invariant  (♦) 
implies  that  (x,y)  is  saturated  after  the  pivot  and  hence  iiuulmissible.  In  this  case  also  it  cannot 
later  be  the  entering  arc  of  a  pivot.  Thus  there  can  only  be  one  pivot  per  arc,  for  a  total  of  at  most 
m.  Once  an  admissible  arc  (x.y)  becomes  a  tree  arc  with  y  the  parent  of  x  and  y  unmarked,  it 
remains  in  this  state  until  it  becomes  the  leaving  arc  of  a  pivot.  It  follows  that  the  total  number  of 
iterations  of  Step  lb  is  at  most  n  +  3m.  □ 

It  is  straightforward  to  implement  Steps  la-lc  so  that  the  total  time  for  Step  1  is  0(nm).  To 
do  this  we  maintain  the  tree  T  using  a  set  of  parent  pointers,  one  for  each  node.  Then  each  itera¬ 


tion  of  Step  lb  or  Ic  takes  0(n)  time. 
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We  can  reduce  the  amortized  time  for  pivoting  and  rerooting  the  tree  by  using  a  dynamic 
tree  data  structure.  The  data  structure  represents  a  collection  of  vertex-disjoint  rooted  tress,  each 
edge  {v.H'}  of  which  has  two  associated  real  values,  g(v,w)  and  g(H',v),  and  each  vertex  of  which 
has  a  mark  bit.  We  denote  by  parentiy)  the  parent  of  vertex  v  in  its  tree;  if  v  is  a  tree  root, 
parent(y)  =  null.  We  adopt  the  convention  that  every  tree  vertex  is  both  an  ancestor  and  a  des¬ 
cendant  of  itself.  The  data  structure  supports  the  following  ten  operations: 

make-tree(y)\  Make  vertex  v  into  a  one-vertex  tree,  with  v  unmarked.  Vertex  v  must  be  in  no 
other  tree. 


find-parentiv):  Return  the  parent  of  vertex  v,  or  null  if  v  is  a  tree  root. 

firul-value(v):  Compute  and  return  g(yparent(v));  if  v  is  a  tree  root,  return  infinity. 

find-min(y):  Find  and  return  an  ancestor  w  of  vertex  v  such  that  g(w,  parendyv))  is  minimum;  if  v 
is  a  tree  root,  return  v. 

find-vertexiy):  Find  and  return  the  unmarked  ancestor  of  vertex  v,  if  any,  nearest  the  root. 

change-markiv,b):  If  bit  6  is  1,  make  vertex  v  marked;  if  b  is  0,  make  v  unmarked. 

change-valueiy,  5):  Add  real  number  5  to  g(yv parent(yv))  and  subtract  5  from  g  {parent{w),w) 
for  every  nonroot  ancestor  w  of  v. 


link(y,w,OL,^):  Combine  the  trees  containing  vertices  v  and  vv  by  making  w  the  parent  of  v. 
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Define  g(v,w)  =  a  and  g(w,v)  =  p.  Vertices  v  and  w  must  be  in  different  trees,  and  v  must 
be  a  tree  root. 

cut(v):  Break  the  tree  containing  vertex  v  into  two  by  deleting  the  edge  joining  v  and  its  parent 
Vertex  v  must  not  be  a  root. 

evert(v):  Reroot  the  tree  containing  vertex  v  by  making  v  the  root. 

We  represent  the  spanning  tree  7  as  a  dynamic  tree.  A  tree  edge  {v.w)  has 
g(v,H')  =  u/(v,w)  and  gfw, v)  =  u/w.v).  (Note  that  we  can  compute  one  of  these  values  from  the 
other,  knowing  u(v,w)  and  u(w,v).)  Initialization  of  7  as  a  dynamic  tree  requires  n  make-tree 
operations  and  n-1  link  operations.  This  initialization  need  be  done  only  once,  at  the  beginning 
of  the  relaxed  network  simplex  algorithm. 

We  perform  Steps  la-lc  using  the  dynamic  tree  operations,  as  follows.  Step  la  requires  n 
change-mark  operations.  There  is  no  need  to  reroot  7;  we  merely  select  as  r  the  current  root  of  7. 

Step  lb  requires  the  ability  to  select  an  admissible  arc  (r,v)  with  v  unmarked,  if  any.  Since 
during  Steps  lb  and  Ic  no  inadmissible  are  becomes  admissible  arxl  no  marked  vertex  becomes 
unmarked,  we  can  perform  this  selection  by  maintaining,  for  each  vertex  x,  a  pointer  into  a  list  of 
incident  arcs  (x,y).  To  find  (r.v),  we  move  the  pointer  for  r  along  its  list  of  incident  arcs,  until 
finding  an  admissible  arc  or  reaching  the  end  of  the  list  Tne  total  time  for  all  such  scans  during 


all  iterations  of  Step  Ib  is  OQn). 
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If  we  do  not  find  an  admissible  arc  (r,v),  we  complete  Step  lb  by  performing  a  change- 
mark  operatioa  If  we  find  an  admissible  arc  (r.v)  that  is  a  tree  arc,  we  complete  Step  lb  by  by 
performing  evertQv).  The  third  and  most  complicated  case  occurs  if  we  find  an  admissible  non- 
tree  arc  (r,v).  Then  we  must  pivot  on  (r.v).  We  compute  x  =  find-min(v)  and  5  =  mm{find- 
valueix),  ujir,v)].  We  increase  the  flow  on  r(r.v)  by  performing  change-valueiv,S)  and  sub¬ 
tracting  6  from  Uf(r,v).  If  ujir,v)  is  now  zero,  we  take  (r.v)  to  be  the  leaving  arc  of  the  pivot; 
Step  lb  is  complete.  If  w/r.v)  is  not  zero,  we  take  (x.  parent  (x))  to  be  the  leaving  arc.  We 
modify  T  appropriately  by  performing  c«/(x)  to  delete  {x,parent(x)},  followed  by 
link  (r.v,  Uf(r,v),  Uf(v,r))  to  add  (r.v).  This  cut  and  link  have  the  side  effect  of  rerooting  the 
tree  at  x,  completing  Step  lb. 

We  perform  Step  Ic  by  scanning  a  list  of  all  the  vertices,  looking  for  an  unmarked  one.  The 
total  time  for  all  such  scans  during  all  iterations  of  Step  Ic  is  0(n).  If  we  find  an  unmarked  ver¬ 
tex  V,  we  letx  =  find-vertex(y),  and  we  complete  Step  Ic  by  performing  evert(x). 

This  method  requires  0(1)  dynamic  tree  operations  for  each  iteration  of  Step  lb  or  Ic,  and 
the  total  time  for  Step  1  is  0(m  logn)  if  the  amortized  time  per  dynamic  tree  operation  is 
O(logn).  The  flows  on  the  tree  arcs  can  be  recovered  at  the  end  of  the  minimum-cost  circulation 
computation  by  performing  n-1  find-value  and  n-1  find-parent  operations. 

We  shall  sketch  an  implementation  of  dynamic  trees  such  that  each  operation  indeed  takes 
OOogn)  amortized  time.  The  implementation  is  that  of  Sleator  and  Taijan  [37],  extended  to 
allow  two  real  values  per  edge  and  a  mark  bit  per  node  instead  of  one  real  value  as  in  [37],  and 
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extended  to  use  reversal  bits  to  handle  evert  operations  as  in  [36].  We  assume  some  familiarity 
with  [37];  we  shall  merely  outline  the  data  structure,  highlighting  the  changes  needed  for  our 
a^lication. 

We  represent  a  dynamic  tree  D  by  a  rooted  virtual  tree  V,  which  contains  the  same  vertices 
as  £>  but  has  different  structure.  Each  vertex  of  V  has  a  left  child  and  a  right  child,  either  or  both 
of  which  can  b  missing,  and  zero  or  mote  middle  children.  We  call  an  edge  of  V  solid  if  it  joins  a 
left  or  right  child  to  its  parent  and  dashed  if  it  joins  a  middle  child  to  its  patent.  Thus  V  corrsists 
of  a  hierarchy  of  binary  trees,  its  solid  subtrees,  interconnected  by  dashed  edges.  The  relation¬ 
ship  between  D  and  V  is  that  the  parent  in  D  of  a  vertex  x  is  the  symmetric-order  successor  of  x  in 
the  solid  subtree  containing  x  in  V,  unless  x  is  last  in  its  solid  subtree,  in  which  case  its  parent  in 
D  is  the  parent  in  V  of  the  root  of  its  solid  subtree.  (See  Figure  1 .)  That  is,  each  solid  subtree  in 
V  corresponds  to  a  path  in  D  from  some  vertex  upward  toward  the  root,  with  symmetric  order  in 
the  solid  subtree  corresponding  to  the  order  along  the  path  from  deepest  to  shallowest  vertex.  We 
say  that  a  vertex  x  is  a  solid  descendant  of  a  vertex  y  in  V',  and  y  is  a  solid  ancestor  of  x,  if  x  is  a 
descendant  of  y  and  the  path  from  x  to  y  consists  of  solid  edges. 

[Figure  1) 

We  represent  the  structure  of  V  by  storing  with  each  vertex  x  pointers  to  its  patent  in  V, 
denoted  by  V-parent(x),  and  pointers  to  its  left  and  right  children,  denoted  by  l^t(x)  and  rightix). 
These  pointers  allow  us  not  only  to  move  up  and  down  through  V  but  also  to  test  in  0(1)  time 
whether  a  vertex  x  is  the  root  of  a  solid  subtree;  such  is  the  case  if  and  only  if  l^V-parentix))  *x 
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and  right(y-parem(x))  ^  X,  We  store  two  bits  with  each  vertex  to  encode  maridng  infonnation. 
For  a  node  x.  bit  mark(x)  is  1  if  x  is  marked;  tot  aU-mark(x)  is  1  if  all  solid  descendants  of  x  are 
maiked. 

We  store  additional  infonnation  with  vertices  to  encode  edge  values.  This  iriformation 
differs  slightly  from  that  used  in  [37]  because  the  values  are  associated  with  edges  of  D  rather 
than  vertices.  We  associate  an  ordered  pair  (x(x).  t(x))  with  each  non/oot  vertex  x  of  V.  as  fol¬ 
lows:  if  X  is  a  left  or  middle  child,  t(x)  is  V-parent(x)  and  s(x)  is  the  last  solid  descendant  of  x  in 
symmetric  order,  if  x  is  a  right  child,  s(x)  is  V-pareroOc)  and  /(x)  is  the  first  solid  descendant  of  x 
in  symmetric  order.  The  correspondence  x  <->  (s(x),r(x)}  is  1-1  between  the  nonroot  vertices 
of  V  and  the  edges  of  D.  We  associate  with  each  nonroot  vertex  x  of  V  four  values: 

g(x)  =  g(j(x),r(x)); 

mg(x)  =  min  {g(y)  |y  is  a  solid  descendant  ofx}; 

h{x)  =  g{t{x),s{x))\ 

mh{x)  =  min  {/i(y)  |  y  is  a  solid  descendant  of  x} 
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We  store  these  values  in  difference  form.  Specihcally,  we  store  the  following  four  values 
with  each  iwnroot  vertex  x: 

(g(x)  if  X  is  a  middle  child; 
g(x)  -  giV-parentiX))  otherwise: 

Amg(x)  =  Ag(x)  -  mg(x); 

(h{x)  if  X  is  a  middle  diild; 

A(x)  -  h(V-parenl(x))  otherwise; 

6mh(x)  =  AA(x)  -  mh(x). 

We  call  the  values  g(x),  Ag(x),  mg(x),  A  mg(x)  the  g-values  of  x  and  define  the  term  h-values 
analogously. 

Finally,  we  store  with  each  vertex  x  a  reversal  bit  rev(x),  whose  interpretation  is  a  follows. 
Let  jrev(x)  be  the  mod-two  sum  of  the  reversal  Wts  of  all  solid  ancestors  of  x.  If  srev(x)  is  1,  the 
true  meaning  of  the  left  and  right  child  pointers  of  x  is  reversed,  as  is  the  meaning  of  the  fields 
holding  the  g-values  and  A-values.  That  is,  lite  left  pointer  points  to  the  right  child  and  vice- 
versa;  the  g-fields  hold  the  A-values  and  vice-versa. 

We  perform  the  dynamic  tree  operations  with  the  aid  of  two  0(l)-iime  restructuring  primi¬ 
tives  on  virtual  trees.  The  first  is  rotation,  in  which  two  vertices  x  and  y  joined  by  a  solid  edge 
are  inteidtanged  while  preserving  symmetric  order.  If  x  is  the  left  (right)  child  of  y  before  the 
rotation,  y  becomes  the  right  Ocft)  child  of  x  after  the  rotation,  and  the  tight  Ocft)  child  of  x 
becomes  the  left  (right)  child  of  y;  all  other  vertices  retain  the  same  parents.  (See  Figure  2.)  It  is 
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straightforward  to  verify  that  all  the  helds  of  all  the  vertices  can  be  updated  in  0(1)  time;  the  only 
subtle  point  is  that  the  ordered  pairs  associated  with  some  of  the  nodes  change.  Specifically,  if 
before  the  rotation  x  is  the  left  child  of  y,  then  the  pair  originally  associated  with  x  becomes  asso¬ 
ciated  with  the  original  right  child  of  x,  the  pair  originally  associated  with  the  original  right  child 
of  X  becomes  associated  with  y,  and  the  pair  originally  associated  with  y  becomes  associated  with 
X.  Nonetheless,  all  the  required  updates  of  values  are  in  vertices  within  two  steps  of  x,  as  are  all 
the  values  needed  for  the  updates. 


[Figure  2] 

The  second  restructuring  primitive  is  splicing,  in  which  a  vertex  y  that  is  the  root  of  a  solid 
subtree  has  its  left  child,  if  any,  changed  to  a  middle  child,  and  possibly  has  a  middle  child 
changed  to  its  left  child.  (See  Figure  3.)  A  splice  requires  less  work  than  a  rotation,  primarily 
because  none  of  the  ordered  pairs  associated  with  the  vertices  change.  This  primitive  also  has  an 
0(1)  time  bound. 


[Figure  3] 

Out  of  rotations  and  splices  we  buUd  the  main  restructuring  operation  on  a  virtual  tree  V, 
called  splaying.  A  splay  at  a  vertex  x  consists  of  a  specific  sequence  of  rotations  and  splices  done 
along  the  path  from  x  to  the  root  of  V.  The  effect  of  the  splaying  is  to  restructure  V,  making  x  the 
root.  The  time  required  for  s[4aying  is  proportional  to  the  original  depth  of  r,  the  amortized  time 


is  OQogn),  if  n  is  the  number  of  vertices  in  V.  The  paper  [37]  gives  the  details  of  how  splaying  is 
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performed  and  the  analysis  leading  to  the  OOogn)  amortized  time  bound. 

We  can  perforai  each  of  the  dynamic  tree  operations  using  at  most  two  splayings  and  0(1) 
additional  tree  restructuring.  We  shall  describe  the  implementation  of  two  of  the  operations;  the 
implementation  of  the  others  is  similar.  To  perform  everi(v),  we  splay  at  v,  perform  a  splice  to 
make  the  left  child  of  v  (if  any)  a  middle  child,  and  flip  the  bit  rev(v).  We  perform  )ind-min(v)  as 
follows.  First,  we  splay  at  v.  Now  the  path  from  v  to  the  root  in  its  dynamic  tree  D  is  represented 
by  the  right  subtree  of  v  in  its  virtual  tree  V.  If  is  the  right  child  of  v,  the  value  OL  =  mg  (x)  is  the 
minimum  of  g(w,  parent(yv))  for  w  an  ancestor  of  v  in  £).  We  walk  down  from  x  through  solid 
descendants,  using  g-values  to  guide  the  search,  until  reaching  a  vertex  y  such  that  g(y)  =  a.  If  y 
is  a  right  child,  our  search  is  complete;  the  vertex  to  be  returned  is  w  =  V-parent(y).  If,  on  the 
other  hand,  y  is  a  left  child,  we  continue  down  the  tree  through  right  children  untD  reaching  a  ver¬ 
tex  w  with  no  right  children.  In  cither  case  we  complete  the  find-min  operation  by  splaying  at  w 
and  returning  w.  (The  splay  pays  for  the  search  to  reach  w.) 

The  dynamic  tree  data  structure  gives  us  an  (7(mlogn)-time  implementation  of  Step  1  of  the 
relaxed  network  simplex  algorithm,  by  Lemma  5.5.  This  result,  combined  with  Theorem  5.3, 
gives  us  our  final  result: 

Theorem  5.6.  The  relaxed  network  simplex  algorithm,  implemented  with  dynamic  trees,  runs  in 
0(nm\ogn  ■  min  (log(nC),  mlogn))time. 
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6.  Remarks 

There  are  two  major  open  questions  related  to  our  work,  one  theoretical,  one  practical.  On 
the  theoretical  side,  we  have  not  resolved  the  question  of  whether  the  networic  simplex  algorithm 
(with  only  cost-decreasing  pivots)  has  a  version  requiring  only  polynomially  many  pivots.  Our 
subexponential  upper  bound  suggests  the  possibility  of  a  nonpolynomial  but  subexponential 
lower  bound.  We  are  not  willing  to  offer  a  conjecture  as  to  whether  the  actual  bound  is  polyno¬ 
mial  or  not,  but  we  believe  that  the  most  fruitAil  way  to  approach  the  problem  is  to  spend  as 
much  time  in  trying  to  construct  a  nonpolynomial  class  of  examples  as  in  trying  to  obtain  a  poly¬ 
nomial  bound.  Whatever  the  answer,  the  results  of  Section  4  show  that  the  key  to  solving  the 
problem  lies  in  obtaining  a  tight  bound  on  the  number  of  pivots  needed  for  cycle  blocking. 

On  tl»e  practical  side,  one  may  ask  whether  any  of  our  results  or  ideas  can  contribute  to 
obtaining  faster  computer  programs  implementing  the  network  simplex  algorithm  or  other  algo- 

t 

lithms.  The  current  computational  wisdom  is  that  the  network  simplex  algorithm  performs  very 
well  in  practice  even  if  no  special  effort  is  made  to  choose  pivots  wisely  [23].  When  advances  in 
computing  technology  make  it  possible  to  attack  problems  on  very  large  networks,  however, 
some  of  our  ideas  may  become  useful. 

In  particular,  the  dynamic  tree  data  structure  can  be  extended  so  that  it  can  be  used  to 
represent  the  spanning  tree  in  the  standard  implementation  of  the  network  simplex  algorithm.  It  is 
necessarily  merely  to  add  values  to  the  vertices  to  represent  arc  costs,  and  to  add  an  additional 
operation  that  sums  arc  costs  along  a  tree  path.  Then  the  amortized  time  to  compute  the  reduced 
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cost  of  an  arc  or  to  perform  a  pivot  is  OOogn).  Standard  implementations  take  0(n)  time  per 
pivot  but  0(1)  time  to  compute  the  reduced  cost  of  r»n  arc.  With  storage  of  one  additional  value 
per  dynamic  tree  vertex,  the  time  to  compute  the  reduced  cost  of  i  arcs,  with  no  intervening 
pivots,  can  be  reduced  to  O  (min{klogn,  k  +  n)).  We  leave  filling  in  the  details  of  this  result  as 
an  exercise. 
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Figure  2.  A  rouUon  in  a  virtual  tree.  Triangles  denote  subtrees. 


